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Relativistic hydrodynamics of an isentropic fluid in a grav- 
itational field is considered as the particular example from the 
family of Lagrangian hydrodynamic-type systems which pos- 
sess an infinite set of integrals of motion due to the symme- 
try of Lagrangian with respect to relabeling of fluid particle 
labels. Flows with fixed topology of the vorticity are inves- 
tigated in quasi-static regime, when deviations of the space- 
time metric and the density of fluid from the corresponding 
equilibrium conflguration are negligibly small. On the base of 
the variational principle for frozen-in vortex lines dynamics, 
the equation of motion for a thin relativistic vortex filament 
is derived in the local induction approximation. 

PACS numbers: 47.75. -ff, 04.20.Fy, 47.15.Ki, 47.32. Cc 

1. In the given work inviscid flows of a relativistic isen- 
tropic fluid are investigated in frame of the general rel- 
ativity theory i.e. in significantly curved space-time 
parameterized by some coordinate system (t, r), with the 
metric tensor gik{t, r). 

For any hydrodynamic system it is possible to separate 
vortex and sound degrees of freedom, with the vortices 
being characterized as the " soft" , while the sound as the 
" hard" degrees of freedom. If there is a static equilibrium 
state in the system with the density of fluid p = po{y) 
and with the velocity v = 0, then the dynamical regime is 
possible with the hard degrees of freedom being excited 
weakly and with dynamics being entirely described by 
slow flows on the background of practically unchanged 
density po(r). A typical velocity in such flows should 
be much smaller then the sound speed (and also smaller 
then the speed of light, as in the case). The general the- 
ory of such slow vortical flows in spatially inhomogeneous 
systems has been developed recently by the present au- 
thor in his work j^]. In particular, the Hamiltonian equa- 
tions of motion have been obtained for an arbitrary given 
topology of the canonical vorticity field, and the general 
form of the Lagrangian for frozen-in vortex lines has been 
established. 

The goal of this Brief Report is to apply the formalism 
developed to the model of relativistic hydrodynamics. As 
the equilibrium state, a spherically symmetric static dis- 
tribution of a fiuid will be considered which is the reason 
for space-time curvature. In the expression for the ac- 
tion of relativistic hydrodynamics the only arbitrary 
dependence s{n) appears which connects the density e 
of the total energy of fiuid measured in the locally co- 
moving reference frame with the density n of number of 
conserved particles The scalar n is equal to the ab- 



solute value of the current 4-vector 

Therefore the static metric of the space-time, as well as 

the equilibrium distribution po (r) of the matter in space, 

after the equation of state e{n) has been fixed, depend 

on single parameter - the total amount of fiuid in the 

system. 

As the result, the equation of motion of a slender vor- 
tex filament in relativistic fiuid will be derived in the 
local induction approximation, and it will become clear 
how the the dynamics of frozen-in vorticity in curved 
space-time differs from the non-relativistic hydrodynam- 
ics, even in the case of small velocities. 

2. Before starting the consideration of the concrete 
dynamical system, it is necessary to make some prelim- 
inary remarks. The ideal relativistic hydrodynamics is 
a member of the large family of the hydrodynamic-type 
Lagrangian models. This family is remarkable because 
of many interesting properties, first of all due to pres- 
ence of infinite number of specific integrals of motion. 
Therefore it is reasonable in the beginning to consider 
the relativistic hydrodynamics from general positions of 
the canonical formalism §,[!). 

As known, the entire Lagrangian description of the fiow 
of some continuous medium can be given by the 3D map- 
ping r = x(i, a), which indicates the space coordinates of 
each medium point labeled by a label a = (oi, 02, 03), at 
an arbitrary moment in time t. The labeling a can be 
chosen in such a manner that the amount of matter in 
a small volume (Pa in the label space is simply equal to 
this volume. Inasmuch as the dynamical system is sup- 
posed to be conservative, the equations of motion for the 
mapping x(i, a) (as well as for the gravitational field, as 
in our case) follow from a variational principle 



SS = sJ /:{x(i,a),5,fc(t, r)}dt^O, 



(1) 



where the Lagrangian £ is a functional of x(i, a), gik{t, r) 
and their derivatives. The explicit expression for C will 
be given later, but now let us note the very important cir- 
cumstance related to the fluidity property of the medium 
under consideration. The fluidity is manifested in the fact 
that the Lagrangian actually contains the dependence on 
x(t, a) only through two Eulerian characteristics of the 
flow, namely through the field of relative density p{t, r) 
and the velocity field v(t, r), i.e. £ = C{p, v, gik}, with 



P(i,r) 



1 



det||9x/5a|| 



(2) 



1 



v{t,r) = x(i,a)|. 



'(*,r) 



(3) 



It follows from these definitions that the dynamics of the 
"density" p{t, r) obeys the continuity equation in its stan- 
dard form 



Pt + V(/9V) 0. 



(4) 



The vanishing condition for variation of the action S = 
J C{p,v, gik}dt when the mapping x(^, a) is varied by 
(5x(i, a), with taking into account the obvious relations 
between the variations 



Sp{r) = -V(p(r) • 5x(a(r))), 



(5v(r) = 5x(a(r)) - ((5x(a(r)) • V)v, 



(5) 



(6) 



can be expressed in Eulerian representation as follows 
(the generalized Euler equation [Q) 



V • V) 



Together with the conditions 6S/Sgik — 0, the equations 
and (Q) determine completely the evolution of hy- 
drodynamic system. It is very important that in all such 
systems an infinite number of conservation laws exists 
[ p2| . The reason of their existence is that the Lagrangian 
£{x(t, a), r)} admits the infinite-parametric sym- 
metry group ~ it assumes the same value on any two 
mappings Xi(i, a) and X2(t, a), if they differ one from 
another only by some relabeling of the labels 

X2(t,a) = xi(t,a*(a)), det||aa7a|| = 1. (8) 

Obviously, such mappings create the same " density" and 
velocity fields. According to the Noether's theorem 

every one-parametric sub-group of the relabeling 
group a* (a) with imit Jacobian corresponds to an inte- 
gral of motion. There are several classifications of these 
conservation laws. For instance, one can postulate that 
the circulation of the canonical momentum p(t, r) 



P(i,r) 



SC 



1 ( 5C{p,w,g,k} 



(5i(a(t,r)) p V 



(9) 



along an arbitrary frozen-in closed contour ^{t) does not 
depend on time (the generalized theorem of Kelvin) 



® (p • dr) — const. 
J -id) 



hit) 

We arrive at a different formulation when consider the 
solenoidal field of the vorticity fl{t, r) 



r2(t,r) = curlp(t,r). 



(10) 



From the equation (|^) the equation of frozenness for 
fl{t, r) follows 



n^ = curl[v X ft]. 
The formal solution of this equation is 
(Oo(a)Va)x(t,a) 



(11) 



0(t,r 



det\\dx{t,a)/da\ 
J (5(r-x(i,a))(f2o(a)Va)x(i,a)da 



(12) 



with the solenoidal field Oo(a) independent on time be- 
ing exactly an integral of motion (the Cauchy invariant 
|l3|). The formula (^2|) displays that the lines of the 
initial field JIq are transported by the fiow, retaining all 
the topological characteristics. Such a property of vortex 
lines is known as their frozenness. 

3. Let's turn our attention to the relativistic hydro- 
dynamics and consider the corresponding expression for 
the action 5 [§, (2) 



S 



1 



g ( £in) + jQ:^^[9ik] ] drdt 



(13) 



Here g = dct\\gik\\ is the determinant of the metric ten- 
sor, R[gik] is the scalar curvature of space-time j^]. In 
this formula, however, the scalar n should be expressed 
through the dynamical variables {p,v, gik}. It is easy 
to see, comparing the continuity equation (^) with the 
equation of fiuid conservation |fi| 



1 d 



-gn- 



dx' 



= 0, 



(14) 



that p 



\f--g dx^ \ ds 
-gn{dt/ds). From this we have the relation 
P 



'9 



goo + goaV^ 



■ gapv^vl^, 



(15) 



which should be substituted into the expression (|l|) and 
after that the equations of motion can be obtained in 
accordance with equations (|^), (^ and SS/Sgik = 0. In 
such form, however, the action contains an arbitrariness 
related to the possibility of coordinate changing. There- 
fore four gauge conditions can be imposed on the metric 
tensor Even more, inasmuch as we are interested in 
a static central-symmetric gravitational field, the corre- 
sponding metric can be chosen as follows 

^^g2ytat ^ ^(^)^^2 _ B{r){dx^ + dy^ + dz^), (16) 



where r = ^y x^ + y'^ + z^. The expression for the scalar 
curvature in this metric takes the form 

1 



- B'A'- 



+AB{A'B' + 2B{2A'/r + A")) 
A'^{-3B'^ + AB{2B'/r + B")) 



(17) 



2 



Therefore the static state is determined by an extremum 
of the functional 



= - / \/AB3 



PO 
/B3 



167r 



r^drdt. 



(18) 



The variations on 5A{r) and on SB{r) should be per- 
formed in the standard way when an extremum is being 
found, while the variation on displacements (5x(a), with 
taking into account the relation and subsequent in- 
tegration, gives the relativistic hydrostatic equation for 
isentropic fluid M 



Po 

/B3 



A = const, 



(19) 



where w{n) = e'{n) is the relativistic enthalpy per one 
conservative particle, while a constant of integration A is 
related to the total amount of matter in the system. 

With fixed static central-symmetric metric (|l^) the La- 
grangian of relativistic hydrodynamics takes the form 



Cr{p,v} 



py/A ~ Bv2 

VaW 



dr. 



(20) 



As far as we examine here only slow flows, we need 
only the first term from the expansion of Cr on powers 
of v"^ 



Po 

/B3 



PqB 



-dr = A 



PqB 
A 



-dr 



Just the functional determines the slow dynamics of 
frozen-in vorticity ft 



CI « curl 



1 SC^ 
Po <5v 



curl I -^^"^ 



4. In order to investigate the motion of frozen-in vor- 
tex lines in the given problem, it is necessary to define 
the Hamiltonian H* js) 



^* " A 



and then express it through the vorticity with help of 
the relation (|l^) and the condition of zero variation of 
density po 



(22) 



The equation of motion for O in quasi-static regime can 
be written after that in the form 



ilf = curl 



curl 



f5ru 

\5n J - poir) 



(23) 



Inasmuch as the vortex lines are frozen into the fluid, it is 
reasonable to consider their shape as the new dynamical 
object 1^ and parameterize the vorticity field as follows 



n{r,t)^J d^i^ S{r~Ii{iy,tt))^idt 



(24) 



Here the label v — {1^1,1^2) G A/" belongs to the 2D mani- 
fold M and singles out a vortex line, while the parameter 
^ determines a point on the line. Obviously, the choice 
of the longitudinal parameter ^ is not unique It 
is easy to understand the meaning of the above formula 
- the frozen-in vorticity field is represented here as the 
continuous distribution of vortex lines. 

It follows from the equations (p3| ) and (|2j) that the 
equation of motion of vortex lines has the Hamiltonian 
form (see 0] for detailed derivation) 



[R,. X Rt]po(R) = 



SH4n{R}} 

(5R 



(25) 



Clearly, this equation does not depend on the choice of 
longitudinal parameterization. It is easy to check by 
direct calculation that the above equation corresponds 
to the variational principle S J C j^dt — with the La- 
grangian Cj^ oi the form |^ 

c^^ j^d^i^ji^ [Kt{v,o X D(R(i.,o)] • Rf(^^,e))rfe 

-KWR}}, (26) 

Here the vector function D(R) is related to the equilib- 
rium density po by the condition [B[ 



(Vr-D(R)) =po(R). 



(27) 



5. Unfortunately, practical use of the new variables 
R(i^, ^, t) is rather difficult in general case because of ne- 
cessity to find the potential component of the canoni- 
cal momentum field p from the system of equations (|lO|) 
and ( p^ ) for substitution of p{r2} into the equation (|2l|). 
However, in analysis of some situations this difficulty 
can be avoided. For instance, the vortex line represen- 
tation allows in the most simple way to study the so 
called local induction approximation (LIA) in dynamics 
of thin vortex filaments. Let's note that in usual Eule- 
rian homogeneous hydrodynamics the LIA gives the inte- 
grable equation which is gauge equivalent to ID nonlinear 
Schroedinger equation [ p^ . Our purpose now is to derive 
the local induction equation for a slender vortex filament 
in a relativistic fluid. Let's suppose that the vorticity is 
present in the system in form of a quasi-one-dimensional 
structure with the circulation F — jj^d^v, with a typi- 
cal width d and with a typical longitudinal scale L ^ d. 
Then we may neglect the dependence of the vortex line 
shape on the label v in the main approximation and con- 
sider the vortex filament as the single curve R(^,i). In 



3 



such conditions, the Hamiltonian is determined by 
the close neighborhood of this curve, so that with the 
logarithmic accuracy is equal to the following expres- 
sion 



(28) 



where the constant A contains the large logarithm 
\n{L/d) > 1 



The density po{r) at r > r* in this case is determined by 
the hydrostatic equation ( p^ with the given A and B. 
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(29) 



With help of the equations (^5|) and (^8|) we can obtain 
the equation of motion for a thin vortex filament in a 
relativistic fluid placed into a central-symmetric static 
gravitational fleld: 



B 



-^d ( 



IRp 



Being solved with respect to the time derivative and 
rewritten in terms of the unit tangent vector t, the unit 
binormal vector b and the curvature k of the line, the 
above equation takes the form 



-(l|)-l^[Vln 



B{R) ) ^ 



Kb, 
(30) 



The condition of its applicability, besides L ^ d, is also 
the inequality 



A_ 
XB 




(31) 



which means the smallness of the maximal velocity in the 
flow in comparison with the speed of light. 

Let us say flnally few words about possible situation 
when the static distribution of the matter in space is 
highly inhomogeneous. Let's suppose that the equation 
of state e(n) results in formation of a dense and massive 
core of relatively small size and extended easy shell 
of a mass being much smaller than the mass of the core 
. Then, as far as we are interested in the dynamics of 
a thin vortex filament in the shell, one can belive that the 
functions A{r) and B{r) in the equation ( |30| ) at r > 
are the same as in the empty space around the mass 
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A{r) = 1 



2r 



(32) 



B(.)=(l + f 



(33) 
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